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ABSTRACT 


The problem of calculating the dynamic stress-strain state of viscoelastic rectangular isotropic plates on a deformed base, including 
those lying freely on the ground medium under the influence of vibration loads, is solved. Several models of the dynamic reaction of 
the base are considered and a qualitative comparison of the results is made. The calculations used the method of Gauss method, 
Muller method and least residuals. 
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1. INTRODUCTION 


Rectangular plates with variable geometric and mechanical parameters under dynamic vibration loads are applied in various 


industries and construction. A rectangular transversely dynamically loaded plate can be supported on a deformable (elastic or 
viscoelastic) base; for example, in the coatings of roads, bridges or runways of air fields. To study the dynamic strength and bearing 
capacity of such structures requires knowledge of their dynamic stress-strain state at vertical loads. The Problem of Flexural 
oscillations of viscoelastic plates on an elastic base is introduced; it is a topical problem of the mathematical theory of viscoelasticity. 
In a closed analytical form, its solution, to simplify the elastic formulation, it is possible to obtain a limited number of boundary value 
problems. An alternative approach to finding an approximate or semi-analytical solution to an elastic problem is to represent the 
solution as a series [1,2,3]. The authors [4] propose, using the variational method for elastic problems, to reduce the resolving 
equations to a system of ordinary differential equations. The disadvantage of these methods is their obvious dependence on the 
methods of setting the boundary conditions and loading patterns. In [4], a finite-difference approach is used for static loading, 
which in turn leads to difficulties in implementing boundary conditions. For high-order differential equations, a large-size template is 
used. All the above reasoning, leads to the need to develop effective methods for solving boundary value problems of the theory of 
plates of the operator on a deformable base (figure 1). 


Figure 1 Plate on elastic base 


2. PROBLEM STATEMENT AND SOLUTION METHODS 


In this paper, the method of collocations and least residuals is used to solve the problem of plate bending numerically. The method 
of finite elements is well established in solving ordinary differential equations and partial differential equations for hydrodynamics 
problems [6]. It is used for the first time to calculate stress-strain state of plates. Consider a rectangular slab on an elastic base. The 
reaction of the elastic base will be considered using a one-parameter model, based on the Winkler hypothesis (hereinafter the 
Winkler model) [5,6,7], and two more complex two-parameter models of Vlasov [4] and Pasternak [8]. Winkler's hypothesis, suggests 
that the reaction of the base is proportional to the deflection of the plate 


p=k,(w-[R@=t)w(e)dr) (1) 
0 


where P) - base reaction, W - plate deflection, k,, - moments coefficient of the bed (the coefficient of proportionality),determined 


experimentally for each type of soil. Despite its simplicity, in many cases the use of this model is sufficient to obtain acceptable 
results from a practical point of view. However, this representation of the soil reaction has a number of drawbacks. For example, 


external loads are distributed on the ground only within the area of the sole of the plate.This position does not correspond to real 
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observations, according to which the soil settles, and therefore is stressed outside the plate.Another disadvantage is the difficulty in 


determining the value of the coefficient k , which depends on the size and shape of the test stamp.A more complex model of soil 


reaction is laid down in two-parameter models 


p=Cy, (w= | R,(t—2)w(2)dzt) -C,,, (Aw = | R,(t-2)Aw(e)d7), 2) 
0 0 


where, A- operator of Laplace, C,,,,C,,,.- soil parameters. Here, in addition to the work of the base on the compression (Winkler 


ml? 
hypothesis), additionally takes into account the work of the base on the shift or slice. In [4], the authors present the base as a 


medium in which there are no longitudinal (along the plane of the resting plate) displacements. Then the coefficients can be 


determined by C,,,,C,,,. the following formulas 


ml? 


E, { E, * 

Cm =z | ORY B.C = dz, 

m Toy Jee are Joe) z 
E Vv 

Ey = - Vo = _— (3) 
l-v, La, 


E, f Ve - moments Young's modulus and Poisson's ratio of the elastic base, Q(z) - the function of the transverse distribution of 


the elastic base, which characterizes the extinction of soil tension with increasing depth H ; 


bz) = er —e!#D) sin (intl) 


where, Y = L3: 


In [8] it is proposed to obtain coefficients from the following considerations. Connects the intensity of the vertical resistance of 


the soil with its sediment, and the second independent coefficient allows you to determine the intensity of the vertical shear force. 
ee fe E,H 


The work also describes the following possible parameters: Cy = ——~ 


:C,.=———. (4) 
1-2." "60 +v,) 


Let's move on to the mathematical formulation of the problem. In a rectangular area O= [0, L] x [0, 1] the boundary value 


problem describing the plate bending taking into account the elastic base reaction is considered (Figure 1) [1,4]. 


2 


D(AAW(x, , X51) — [ R@—r)AAw(x,,x,,7)dr) + pH on = Q(X, %_,t) — p(%,,%,t), 


>= 
where W(X,,X5,f) -the deflection plate; g(X,,X>,1) - external load; p(X,,X,,1) - the reaction of the elastic foundation; 


D a a /42(1—v’)) - cylindrical stiffness;l,,1,,h - length, width, thickness of the plate; £), V - moments Young's 


modulus and the Poisson's ratio of the plate. 


The elastic base reaction is determined for each model from the corresponding formulas (1), (2) with coefficients (3) or (4). By 


P(x, x.) = 0 we obtain the classical equation of plate bending [1]. On the edges of the plate we will use the known boundary 


conditions [1]. For example, when X, = 0 maybe a free edge: 
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2 2 3 3 
CALA ALAN YS Ye 
Ox, Ox; Ox; Ox,Ox; 


=Q!. 


Special attention should be paid to the size Q! . This function can be interpreted as the influence of soil outside the plate on its 
edges [4,8]. Since the Winkler model does not account for this effect, then for it O! =(). For the two-parametric models of the 


edges Q! takes the following form [4] 


2 
Ow 1/( Ow 
Q° =C,| aw+——-—| — 
Ox, 2a\ Ox; 
Similarly, you can write conditions on the other edges of the plate. 


Let's cover Q the area with a rectangular grid with cells uniform in each direction Q, (I =... N ). To determine the solution 


in each cell, we will use the domain decomposition method-the method of iterations on subdomains (alternating Schwartz method), 


in which the subdomain is a cell. In each cell, you enter a local coordinate system associated with the source variables by the 


* * 
following formulas y= (x, xX, )/h, y= (x, xX, )/h, , where 2h, ,2h, - the size of the cell in the direction Xy, 


ok * 
X, respectively; (x, ee - cell center coordinate. In each cell, the approximate solution is represented as a polynomial of the 


fourth degree and to determine the unknown coefficients, we write a local system of linear algebraic equations . This system 


includes collocation equations 


he d*w Ow ko we | h2 Otw'(r) Otw* (tr) hh? O4wi (rc 
Dy( q “whl a : 4 frie 7)( : a di = 5 x yar) 
h, Oy, Oy, Oy, hy Oy, = h, oy, Oy, Oy, hy Oy, 
O°w* O° w* 0’w 
-C,, (h; Lh . )+hoheC,, wt + 0H =h?h’¢q, 
2, ay? 1 ay: 2h Cin Pp ar 2h | 


k 
where W, - shuffle the point touxu k - the body, 1 - the outer normal to the border @. . The boundary condition has the form in 


case of pinching 


The local system of linear algebraic equations consists of 9 integro-differential collocation equations written at the inner points 
of the cell. Also, at each cell boundary, depending on whether this boundary is adjacent to the boundary of the source region, three 
matching conditions or three boundary conditions are written. The resulting equations will be overridden. Its solution will be 


understood in the sense of least squares. 


3. NUMERICAL RESULTS 


Consider a rectangular plate on an elastic base under the action of a uniform dynamic load g = Q,e [9,10] (figure 2 & 3). 
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Figure 2 The amplitude Shape of a uniformly loaded plate whose two edges are pinched for Winkler (a), Vlasov (6), Pasternak (B) 
models. 


Figure 3 Section of the plate deflection amplitude at, two edges of which are pinched, for Winkler models (solid), Vlasov (dashed), 
Pasternak (dashed). 


Two adjacent sides of the plate are pinched, the other two are free. In the experiment, the calculations for three models of the 


base (Fig. 2, 3) for parameters iP = 2s = 20 M, i= 0.1m, H=2 M, Eo=200INa, V = 0.28, E, =O 4 Ma, Ve = 0.4, 


k =0.31rna/m, Qo=1 MMa. It can be seen from the figures that for two-parameter models, taking into account the function on 
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the free edge leads to its lifting, which from the point of view of real experience is more logical than for the case of the Winkler 


model, when the free edge is deformed without bending. 
Consider a square plate, free-lying on the ground, modeling, for example, the Foundation of the bridge support. The plate is 


under the action of a uniform dynamic (harmonic) load applied to the region [2,8] x [2,8] (rice. 4, 5), L, =I, = 10m, 
h=0.1m, H =2m, E=200rnav =0.28, E, = 0.4rna, Ve =0.4, k =0.3rna/m, Qo=1 MMa. In this case, the 


deflection of the plates does not depend qualitatively on the choice of the base reaction model, since the values of deflections on 
the contour are small (figure 4 & 5). 


Figure 4 The amplitude Shape of the free-lying plate on the ground under a special load for the Winkler (a), Vlasov (6), Pasternak (B) 
models. 


Figure 5 Section of the deflection amplitude of the free-lying plate on the ground for Winkler models (solid), Vlasov (dashed), 
Pasternak (dashed). 
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4. CALCULATIONS 


Thus, an algorithm for solving dynamic and quasistatic problems is given in the paper. It was established that the models of the 


Winkler and Pasternak base give almost the same results, their difference is up to 5-7%, and the difference in the models of the 
Vlasov base is up to 30%. This can be seen from the figure. 5. It was also found that the deflection of the plates does not 
qualitatively depend on the choice of the base reaction model, since the deflection values on the circuit for all the models 
considered differ up to 10%. 
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